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Abstract 

A prototype for a sampling calorimeter made out of cerium fluoride crystals interleaved with tungsten 
plates, and read out by wavelength-shifting fibres, has been exposed to beams of electrons with ener¬ 
gies between 20 and 150 GeV, produced by the CERN Super Proton Synchrotron accelerator complex. 
The performance of the prototype is presented and compared to that of a Geant 4 simulation of the 
apparatus. Particular emphasis is given to the response uniformity across the channel front face, and 
to the prototype’s energy resolution. 
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1 Introduction 

Future hadron colliders are expected to have high instantaneous luminosity (> 10^^ cm“^ s“^) and will operate 
for long data-taking periods, delivering to experiments large quantities of data (> 1 ab~^). This high-radiation 
environment sets stringent requirements on the radiation hardness of all detector parts, and on the components 
placed in the forward region in particular. 

A sampling calorimeter made of cerium fluoride (CeFa) scintillating crystals, interleaved with tungsten plates, 
constitutes a viable option for forward electromagnetic calorimetry. Cerium fluoride spontaneously recovers, at 
room temperature, from hadronic damage m, and its stoichiometry can be tuned in order to make it extremely 
resistant to ionizing radiation m. These features, together with the fact that its scintillation light has a spectrum 
which is suitable for wavelength shifting (the scintillation peak is around 340 nm) and is fast enough for the high- 
repetition frequencies of modem colliders (it has a decay time of 30 ns), make cerium fluoride a very attractive 
active material. 

Stemming from the design of a previously tested prototype a, a new design has been built, consisting of a single 
tower of 15 layers of alternating CeFa crystals and tungsten absorber plates. This article describes the results of 
irradiating this prototype with electrons of energies between 20 and 150 GeV, produced by the CERN Super Proton 
Synchrotron (SPS) accelerator complex. 


2 The W-CeFa Prototype 

A picture of the prototype during assembly can be seen in Figure [T] (left). It is made of 15 layers of 10 mm-thick 
CeFa crystals, interleaved with 3.1 mm-thick tungsten absorber plates, for a total length of 25 Xq. The Barium- 
doped CeFa crystals have been produced by Tokuyama a, and their faces have been polished for total internal 
reflection. The module transverse dimension of 24 x 24 mm^ matches the effective Moliere radius of the sampling 
compound, which is found to be 23.1 mm from a Geant4 Q simulation. A higher granularity can be achieved 
by increasing the absorber fraction, at the cost of lower sampling fraction. 

Four 3 mm-wide chamfers are cut on the short edges, allowing wavelength-shifting (WLS) fibres to run alongside 
the tower, transmitting the scintillation light to the back of the detector. The chamfers have been lapped to a 
roughness of Ra = 0.4 ± 0.1 /rm, so as to allow the light to escape the crystals and be collected by the WLS fibres. 
A 0.2 mm layer of highly-diffusive DuPont™Tyvek® |j^ foil is inserted between each crystal and its neighboring 
absorber plates, in order to maximise the light collection of the fibres. 

In this beam test four 3HF fibres produced by Kuraray 13 have been employed. The fibres have a peak excitation 
of 340 nm, compatible with the emission spectrum of the CeFa, and an emission peak at 540 nm. The fibres have a 
single 0.02 mm layer of cladding made of Polymethyl Methacrylate (PMMA) surrounding the core of Polystyrene. 
One end is coupled with optical couplant to a Hamamatsu m R1450 photomultiplier tube (PMT) at the rear end; 
the other end is aluminized. The quantum efficiency of the PMTs for the light emitted by the fibres has been 
measured to be about 7%. The radiation hardness of these fibres is not sufficient for future applications at the 
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HL-LHC, therefore R&D targeted at developing radiation-hard fibres is currently ongoing. 

For lateral shower containment, the W-CeFa prototype has been surrounded by 24 BGO crystals, to complete a 
5x5 matrix. The BGO crystals have been taken from the disassembly of the electromagnetic calorimeter of the 
L3 experiment operating at LEP, and are 24 cm long and have a front face of 22 x 22 mm^ and a rear face 
of 30 X 30 mm^. The light from each BGO crystal is read out with a Hamamatsu R1450 PMT. A photograph of 
the front of the assembled 5x5 matrix is found in Figure(right). The matrix was put in a light-tight box, and 
kept at a temperature of 18.0 ± 0.1° C. The box was placed on a moving table, which allowed to center the beam 
anywhere on the matrix, with a precision of 0.1 mm. 

The fibre signals were then read out by a CAEN VI742 digitizer Qol, while the BGOs were read out by a CAEN 
V792 ADC. The digitizer gives access to the full pulse-shape of the signal, which is stored in 1024 samples, at 
a frequency of 5 GHz. The width of the readout window allowed a baseline subtraction on a per-event basis. A 
low-pass filter, which was found to reduce high-frequency noise present in the signal and improve the stability of 
the measurement, was applied to the output of the digitizer. The channel response is defined as the full charge- 
integrated signal, after the baseline subtraction. 


3 Experimental Setup and Electron Selection 

The high-energy electrons from the H4 beam at the CERN SPS North Area are tracked before hitting the front face 
of the W-CeEa prototype: a wire-chamber (placed at 12 m from the channel), two scintillating-fibre hodoscopes (at 
6 and 3 m, respectively), and a small hodoscope, just in front of tower front face. 

The wire chamber is made of two planes of 55 cathode wires and 28 anode wires, organized in a grid with an active 
area of 80 x 80 mm^. Each plane provides, respectively, a measurement of the electron position in the x and y 
directions with a nominal resolution better than 200 pm im. 

Each of the two scintillating-fibre hodoscopes is composed of two layers of 64 plastic fibres of 0.5 mm diameter, 
oriented in the x and y directions, respectively. The signals from the fibres are clustered grouping together adjacent 
fibres, up to a maximum of 4 fibres per cluster. The cluster position is defined as the average position of its fibres 
and it is used to estimate the trajectory of the particle before impacting the calorimeter. 

Einally, a small 2x2 mm^ hodoscope is placed in front of the calorimeter front face, and aligned with the center 
of the CeEs module. This hodoscope is made of two pairs of 1-mm fibres, again oriented in the x — y directions. 
Being fixed to the calorimeter box, this small hodoscope is used to align the center of the tower front face with 
respect to the tracking detectors: events are selected in which the particle produces a signal in one of the two fibres 
of the small hodoscope, and the average recorded position is measured in each of the other three tracking devices. 
The average position measured by each tracking device is then used as an alignment offset, and is subtracted 
from subsequent measurements. This is done separately for each beam energy. After the alignment procedure the 
hodoscopes and the wire chamber are aligned within 0.5 mm. 

Eour plastic scintillators of varying size, the smallest of which has a transverse dimension of 1 x 1 cm^, are placed 
on the beamline. The data acquisition is triggered by the coincidence of all four scintillators. 

Beam electrons are found to most likely produce 2-fibre clusters. Electrons hitting upstream material will pro¬ 
duce particle showers before hitting the tower; these events are vetoed by selecting only events which registered 
hodoscope activity compatible with the passage of one electron. It is also additionally required that electrons are 
parallel to the beam line by requiring that the position reconstruction of the wire chamber and the far hodoscope 
are within 4 mm of each other, and that the positions measured in the two hodoscope planes lie within 1.5 mm. The 
reconstructed impact position of the electron on the calorimeter is taken as the average of the positions measured 
by the two hodoscopes. 


4 Experimental Setup Simulation 

A full simulation of the experimental setup has been implemented within the Geant4 framework. The W-CeEa 
channel description has a realistic geometry of the CeEa, tungsten, and the fibres. The BGO crystals are approx¬ 
imated by 22 mm x 22 mm x 24 cm parallelepipeds. The hodoscopes and scintillators along the beam line are 
approximated with blocks of PMMA with the corresponding sizes. 

The simulated response of the channel is taken as the sum of the energy deposited in the active material, and 
each layer is weighted by the relative light yield difference of the 15 crystals, which has been measured in the 
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100 GeV Electron Beam 



Figure 2: Example response distribution, for 100 GeV electrons. The result of a Crystal-Ball fit is superimposed 
with a blue line, and the mean p and gaussian width a of the fitted function are shown. 
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Figure 3: Response of the W-CeFa channel as a function of the electron impact point on the front face for 
50 GeV (left) and 100 GeV (right) electrons. 

laboratory with a radioactive source, and found to have a spread of 9%. The simulation automatically takes into 
account resolution effects due to stochastic shower fluctuations and lateral containment. An additional term to 
include photostatistics, determined by the in situ absolute photoelectron calibration from cosmic muons, is added 
in quadrature. 

In addition to the experimental setup, an additional setup, consisting of a 5 x 5 matrix of W-CeFa towers, has 
been simulated. This allows for a more realistic description of the lateral leakage and for an extrapolation to future 
calorimetric applications. 


5 Results 

The response of the four read-out channels of the prototype is intercalibrated in order to take into account differ¬ 
ences in the gains and quantum efficiencies of the PMTs. This is done by selecting events where electrons impinge 
on a narrow ±3 mm strip passing through the center of the front face, alternatively in the x and y directions, and 
studying the variation of the response of each single read-out channel as a function of the electron impact position 
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Figure 4; Response profile as a function of the beam impact position, projected along the x (left) and y (right) 
directions, for data (solid markers) and the simulation (open markers). 

along that strip. The inverse of the response at the center of the strip is taken as intercalibration factor, and the 
two factors obtained in the x and y directions are averaged. The resulting intercalibration factors are applied as a 
multiplicative, relative, correction. The four calibration factors are found to be within 10% of unity. 

The response of the W-CeFa channel is defined through a fit to the sum of the inter-calibrated signals recorded 
by the four fibres. The adopted fit function is a Crystal Ball uni [El, to model the low-energy response tail. An 
example fit, for 100 GeV electrons, can be seen in Figure]^ 

The resolution of the tracking devices on the beam allows for a precise study of the prototype’s response as a func¬ 
tion of the electron impact point on the front face. This is shown in Figurej^for 50 GeV (left) and 100 GeV (right) 
electrons. The nominal position of the channel’s front face is overlayed as a yellow frame, and the agreement with 
the observed response profile is proof of a successful alignment. Furthermore, the response is found to be quite 
uniform across the channel’s front face. 

Figure|^shows instead the response profile as a function of the beam impact position, projected along the x (left) 
and y (right) directions, respectively for data (solid markers) and the simulation (open markers). The simulation 
does not include ray tracing and light collection effect. Nevertheless the difference between the data and simulation 
is at most 5%, so we can conclude that the effect of light collection in the four fibers is smaller than 5% in this 
configuration. In order to minimize the effect of lateral containment fluctuations, for the resolution studies that 
follow, events impacting the central 6x6 mm^ area of the channel front face have been selected. 

The resolution of the response distribution (jr is defined as the ratio between the mean and the width of the gaussian 
core of the Crystal Ball fit. To extract the energy resolution as of the W-CeFs tower, it is necessary to disentagle 
(i) PMT saturation effects (present at high voltage settings) and (ii) PMT noise (present at low voltage settings). 
Because of the dependence of both effects on the high voltage applied to the PMTs, in situ voltage scans have been 
taken at different beam energies. 

The result of these scans is shown in Figure]^ where cjr is studied as a function of the mean response R, for beam 
energies of 20 GeV (squares), 50 GeV (circles), 100 GeV (stars) and 150 GeV (crosses). The points are fitted with 
the function: 

/ X N 
(Jr[R) = 

where the first term describes PMT noise and is governed by the parameter N which is fitted simultaneously to 
all energy response scans, and the second term is the instrinsic energy resolution, at a given energy. The results 
of the fits are superimposed in Figure]^ The uncertainties of the points shown in Figure [^include a systematic 
uncertainty related to the choice of the functional form chosen to fit the response distributions, which has been 
evaluated separately for each energy, and reaches a maximum value of 20% at the lowest energy point. 

In order to limit the effect of PMT saturation, which manifests itself at high response, the fitting range is limited to 
0.7 X 10® ADC channels. The fit is repeated twice, with an upper limit of 0.5 x 10® and 0.9 x 10® ADC channels, 
respectively, and the maximal variation in cr^; is taken as a systematic uncertainty on the choice of the fitting range. 
It is found to be less than 10% for all energies. 

The resulting energy resolutions of the prototype are shown with solid blue markers in Figure]^ as a function 
of the electron beam energy. The data are compared to the results obtained with the simulation (hollow markers) 
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50 GeV Electron Beam 
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Figure 5; Response resolution gr as a function of the response, for beam energies of 20 GeV (squares), 
50 GeV (circles), 100 GeV (stars) and 150 GeV (crosses). The results of the fits are superimposed. 
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Figure 6: Energy resolution of the W-CeFa single tower as a function of the electron beam energy (solid blue 
markers), compared to the simulation (hollow markers). Resolution of the simulated 5x5 W-CeFs matrix is 
shown with square black markers. 

and the two are found to be in agreement. This allows for an extrapolation to the simulated 5x5 simulated matrix, 
whose results are shown with black square markers, and are fitted with the function; 
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where S is the stochastic term and C is the constant term. We obtain S = 9.54 ± 0.13% and C = 0.33 ± 0.04%. 
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6 Conclusions 


A prototype for a sampling calorimeter made out of cerium fluoride crystals interleaved with tungsten plates, and 
read out by wavelength-shifting fibres, has been exposed to beams of electrons with energies between 20 and 
150 GeV, produced by the CERN Super Proton Synchrotron accelerator complex. The tower response is found to 
be uniform across the channel front face. We find the measured energy resolution to be in good agreement with a 
dedicated Geant4 simulation. According to the same simulation, a 5 x 5 tower matrix would achieve a stochastic 
term of (9.54 ± 0.13 )%/-^/E [GeV]. This constitutes a proof of concept of a technology with potential use in 
future hadron-collider applications. 
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